1 

Performance analysis of an opportunistic relay 
selection protocol for multi-hop networks 

Technical report 

Kostas Stamatiou, Davide Chiarotto, Federico Librino and Michele Zorzi 
Department of Information Engineering, University of Padova 

Abstract 

In this technical report, we analyze the performance of an interference-aware opportunistic relay selection protocol 
for multi-hop line networks which is based on the following simple rule: a node always transmits if it has a packet, 
except when its successive node on the line is transmitting. We derive analytically the saturation throughput and the 
end-to-end delay for two and three hop networks, and present simulation results for higher numbers of hops. In the 
case of three hops, we determine the throughput-optimal relay positions. 

I. Introduction 

Opportunistic routing in multi-hop wireless networks takes advantage of favorable channel conditions in order 
to advance packets over a large distance, thus reducing the end-to-end packet delay. In this paper, we consider 
a line network consisting of a source, a number of relays and a destination, and evaluate the performance of an 
"interference-aware" opportunistic relaying protocol; namely, a node always attempts to transmit a packet to the 
farthest node down the line, except when its successive node is transmitting, in which case it stays silent. The 
rationale of the protocol is simple: if two consecutive nodes transmit, it is unlikely that the transmission of the 
node farther from the destination (FAR) will be successful, due to the strong interference generated from the node 
closer to the destination (CLOSE); therefore, FAR stays silent in order to avoid interfering with the transmission 
of CLOSE. 

In the case of a two-hop system, we derive exact expressions for the saturation throughput and the mean end- 
to-end delay, while in the case of three hops, we obtain an exact expression for the saturation throughput and an 
accurate approximation for the mean end-to-end delay. The analysis takes fully into account the interaction between 
the source and relay queues and is based on a generating function approach employed in early work on packet radio 
networks JT). To the best of our knowledge, these are the first analytical expressions for the delay and throughput 
of tandem queueing networks with opportunistic routing and a realistic underlying physical layer model that takes 
into account fading and interference. 
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We provide numerical and simulation results for a path-loss and Rayleigh fading channel. In particular, for the 
case of three hops, we determine the relay positions that maximize the saturation throughput. Overall, for typical 
values of the average link signal-to-noise ratio (SNR), the throughput gain of the considered protocol with respect 
to an aggressive opportunistic relaying protocol is 10-15%, and even larger with respect to a TDM A protocol. 
Simulation results for four and five hop systems exhibit similar performance gains. 

Early work on tandem queueing networks [2] relied on simplified channel and interference modeling and did 
not consider direct packet transmissions over the distance of multiple hops. Recent work on opportunistic routing 
includes [3|-|5|. Various aspects of line networks have been studied in [6|-[8|, while [9] calculated the end-to-end 
throughput of dynamic relay selection in a random geometric setting. 

II. System model 

We consider a slotted-time system, where slot t e Z is the time interval [t, t + 1), and the slot duration, i.e., 
unity, is equal to the duration of a packet. The system consists of N + 1 nodes, i.e., the source, N — 1 relays and 
the destination. At the end of each slot, a new packet arrives at the end of the source queue with probability A and 
arrivals are independent across slots (other arrival distributions can also be accommodated by the analysis). The 
buffer size at the source is infinite. According to the considered protocol, a node transmits its head-of-line packet 
in slot t, if its successive node does not transmit in that slot (the last relay always transmits since the destination 
acts as a sink). The packet is kept at the farthest receiver that successfully receives the packet, and is discarded 
by all others. If the packet is not successfully received by any receiver, it remains at the head-of-line. Finally, we 
assume that nodes can not transmit and receive simultaneously. 

For analytical purposes, we make the following assumptions: 

• (Al) A packet can cover the distance of at most two hops; 

• (A2) interference from a transmitter more than two hops away from a receiving node is negligible; 

• (A3) the buffer size at the relays is unity. 

(Al) and (A2) are based on the fact that, in terrestial networks, the signal power decreases quickly with distance 
due to the large path-loss exponent. Therefore, a direct three-hop transmission is highly unlikely for typical SNR 
values and the interference from far-away transmitters is close to negligible. These statements are also justified by 
the simulation results of Section [TV] Regarding (A3), a relay buffer size larger than unity is unnecessary for N = 2, 
since, by virtue of the protocol, the only relay will always transmit if it has a packet, thus it can not receive. For 
N > 3, a buffer size larger than unity could enable a relay to receive a packet in the event that its successive relay 
transmits. Nevertheless, in Section [IV] it is demonstrated via simulation that the protocol performance is insensitive 
to the relay buffer size for three, four and five hop systems. 

Let the numbers 1, . . . , N correspond to the source, the first relay,. . . , the (N — l) th relay. In the absence of 
interference, we denote the probability that a transmission of node n succeeds in covering two hops and one hop as 
P20,n and pio,n, where pio, n > P20,n- (we set P20.N — since the last relay cannot perform two-hop transmissions). 
The probability that a packet covers at least one hop is p s _ n = pio,n + (1 — Pio,n)P20,n (where the subscript "s" 
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stands for "success"). We also define n — 1, . . . , N — 2 as the probability of successful reception over a single 

hop, in the presence of interference from transmitter n + 2, i.e., one hop away from the receiver of n. Henceforth, 
we employ the following notation. The complement of x is x, i.e., x = 1 — x; the derivative of the function f x is 



f or ^ 

Jx U1 dx 

|A|. 



f x \ the double derivative of f x is /" or -n^f 



/" ; and the determinant of matrix A is 



III. Analysis 

Let Qi(t), Qi(t), ■ ■ ■ , Qiv(i) denote the number of packets at the source, the 1 st relay,. . . , the (N ~ l) th relay, 
respectively. Since the relay buffer size is unity, Q n (t) € {0, 1} for n = 2, . . . , N, while, for the source, Qi(t) € N. 
In steady state, the probability generating function (pgf) of the vector (Qi(t), . . . , Qjv(t)) is 



9x 1 X2--.X N 



r QiW Qa(t) 



r Q«(t) 



(1) 



where (x\, . . . , xjv) € [0, 1] N is the argument of the pgf. Moreover, let A(t) be a Bernoulli random variable with 
parameter A which represents the arrival (or not) of a new packet at the source at the end of slot t. Then, the pgf 



of A{t) is f Xl = E 



A(t) 



= Xxi + A. 



The mean end-to-end delay is calculated as [2| 



D 



1 N 

IS 



dg Xl ... Xf 



n=l 



dx n 



-x N = l 



The saturation throughput t s is defined as the minimum value of A for which D becomes infinite. 



(2) 



A. Two-hop network (N = 2) 

In steady state, g xlX2 satisfies the equation 



9X! 



h(t+l) Q 2 (t+l) 



Q 2 (*+i) 



. 2 (I(Qi(i) = 0,Q 2 (t) = 0)+I(Qi(t) >0,g 2 (<) = 0) 
+ I (Qi(t) = 0, Q 2 {t) > 0) + I (Q x (t) > o, g 2 (t) > 0))] , 



(3) 



where !(•) is the indicator function. From (0 and ([T), it follows that g XlX2 must satisfy the functional equation 



9x!X 2 fxi 



.9oo + P10P20 — I 1-P10P20 \9xi_o - goo) + l-pio {gx 1 x 2 - gxxO) 

V Xx Xl J \x 2 ' 



(4) 



The first, second and third terms in the brackets corresponds to the following events: both the source and relay are 
empty; only the source is non-empty, thus a packet advances directly to the destination with probability p 2 o, or to 
the relay with probability P10P20, or to neither with probability P10P20', the relay is non-empty (and the source is 
either empty, or remains silent if it is non-empty), thus the packet transmission to the destination succeeds with 
probability p w or fails with probability p w . 

In Proposition Q] we derive the delay and saturation throughput of a symmetrical two-hop network, i.e., p±o,i = 
Pio,2 = Pio- Since only the source can perform a two-hop transmission, for simplicity we write p 2 o,i = P20- 
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Proposition 1 For a symmetrical two-hop network, the mean end-to-end delay is 



1_A 1-212220 , 



ana! the saturation throughput is 

Ps 



1 + P20 

Proof: Since the size of the relay buffer is unity, from ([!}, we have that 



Solving this system of equations with respect to g Xl o, gx x i, we obtain 

3 0= goof ^lihl Z l )+PWP4 1 ~ Xl)f Xl - (P20 + PlOP2QXl){f Xl - 1) ^ 

""" Xl Xl{f Xl - 1) +Pl0Ps(l - X\)f% x - (P20 + (p~10P~20 ~ Plo)xi) fxdfx! ~ 1) 



A \ dxi 



Xi — l 



(6) 



g XlX2 = g Xl o + (Sail - 5no) ^2- (7) 

Substituting in ©, we obtain 

9xix2 f^ 1 = 9oo + {P10P20X2X1 1 +P20X1 1 + P10P20) {9xi0 ~ 9oo) + (P10+P10X2) (g xi i ~ flWo) ■ (8) 
Letting £2 = {0, 1} in © yields 

feO (l - fx! (JP20X1 1 +P10P20 -Pio)) ~ 9x-LlfxxPW = 500./W (l ~ P20X^ ~ PlOP~2o) 

gx 1 ofx 1 Ps{x^[ 1 - l)+g Xl \(f Xl - 1) = soo/xiP^^r 1 - !)• 



= Qo /W (1 - xi) (xi(l - f xi ) +p w xif xi ) 

S sci afi(/an - 1) + PioPs(l - x 1 )fl i - (p 2 o + (P10P20 - Pio)zi) /n(/n - 1) 
Letting x — > 1 in (TlOt and applying de l'Hopital's rule results in 

, x 1 + P~20 

goo = 1 - A . 

Ps 

From (O and (|7), we have 

D = U d 4^ + i-g w ), (ID 



where the first and second terms in the parentheses are the mean queue sizes at the source and relay buffers, 
respectively (the latter is equal to the probability that the buffer is not empty, since the buffer has size unity). From 
<f9t-(fT0b. with the help of de l'Hopital's rule, we determine 510 and . After some algebra, we obtain (|5j. 

x — l 

Eq. (O follows from the definition of t s . ■ 
As seen in (O, the packet arrival rate which saturates the source buffer is given by (pio +P2o(l —Pio))/ (2 — P2o)- 
The expression clearly shows the gain of opportunistic routing (P20 > 0), with respect to a protocol where two-hop 
transmissions are not allowed (P20 = 0), in which case r s = pio/2. 
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B. Three-hop network (N = 3) 

Since pn. n is defined only for n = 1, for simplicity we write p\\ \ — p\\. The system generating function 



satisfies 



fxi 



ffOOO H (9x 1 x 2 x 3 — g Xl x 2 Q — 9x 1 0x 3 + fcOO + <?00x 3 - 5OOOJ H (<?XiX 2 ~ ffxiOO) 

^3 ^2 



(ff^ioo — 5ooo) 



■^3 



"(<feiOx 3 — fcOO ~ 300^3 + 300o) 



(12) 



where 



a X3 = Pl0,3 +^10,3^3 

i 2 i 3 = PlO, 2^20,2^2 + PXO, 2P20,2a;3 + P20,2 

-1 



Cxix 2 x 3 =PlQ,lP20,l +^10,1^20,1^22:1 



KO,!^! 



(piiX2a; 1 1 + pn)(pio,3 + PlO^)- 



(13) 



The different terms on the right hand side of (TTZt can be explained in a fashion similar to ©. The main difference 
between the two equations is the presence of the last term in (fT2l . which captures the event of concurrent 
transmissions from the source to the first relay and the second relay to the destination. In Proposition [2] we 
derive the end-to-end delay and the saturation throughput of a three-hop network. 



Proposition 2 For a three-hop network, the end-to-end delay is 

D 

and the saturation throughput is 
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(16) 
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and all derivatives in M4\) -M5\l are taken with respect to X\. The constant C is defined as 

£ a .9poi _ 2 
5ooo 

Proof: Since the size of the relay buffers is unity, we have 

5xix 2 x 3 = 5xix 2 "I" (5xix 2 l 5xix 2 o) x 3 
5xix 2 x 3 = 9xi0x 3 ~\~ (5xilx 3 9x 1 0x 3 ) x 2- 

Therefore, (TTZt becomes 

9x!x 2 x 3 — fx\ [5000 + a x 3 (9x 1 x 2 i — g Xl x 2 o — g Xl oi + SxiOO + ffooi — ffooo) + b X2X3 (g Xl 10 — 5xiOo) 

~r C-X1X2X3 (ftriOO — ffOOo) + ^xiX2X 3 (SxiOl — 5xi00 — 5001 + 500o)] ■ 

Setting (x2,xs) — {(0, 0), (0, 1), (1, 0), (1, 1)} in (|22T > and x\ = x (to make the notation easier), we obtain 

boiSxio + (dxoi — f x 1 )g x oi + (— boi + c x oi — d x oi)5xoo == ( — C — 1 + c^oi + Cd x oi)gooo 
a-ogxii + (—0,0 — fx 1 + bio)<7xio + ( — fl o + d x io)g x oi + (ao — bio + c x io — dxio^xoo = 

(— 1 — OoC + c a;10 + Cda;lo)5000 

booffxio + d x oog x oi + ( — fx 1 — boo + c x oo — dxoo).9xoo = ( — 1 — a oC + c :roo + C^xoo)ffooo 
(1 - f x )9xii + (-1 + d xl i)g x0 i + (cu - dxii)ffx00 = (-C - 1 + c x \\ + Cd x u)gooo, 
where the constant C is defined in ( f2Qb . Solving (f23b-(f26b over g'xiii^xiOjfl'xOijfl'xOO) we obtain 
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where K x n, K xW , K x0 i, K x00 are defined in (O. 



-j 5x0i — 5ooo- 



xOl 



K 



Kr 



5x00 — 5ooo- 



j:00 



Kx ' 



7 



By the law of total probability, gm = 1. Therefore, taking the limit of g x n for x — > 1 and applying de l'Hopital's 
rule, we have 



ffooo — 



where 



boi dioi — 1 cioi — 1 
boo dioa cioo — 1 
~d'm — Cm 



boi d\oi ~~ 1 c ioi — 1 

bio - oo - 1 cno — oq — 1 
boo ^ioo cioo — 1 



(28) 



(29) 



and 

Kin = a Q 



boi d' in (cioi - 1) - c' m (dioi - 1) 
boo d'm( c ioo — 1) — c'ih^ioo 



+ a C 



b i d' m (cioi - 1) - c' m (dioi - 1) 

i>io — do — 1 rf' m (ciio - a - 1) 



(30) 

From the definitions in ( fT3l , it is straightforward to show that all the determinants in d29li-(|30li are positive. Moreover, 
C > 0, since 5o0 i = P(Qi(i) = Qa(*) = 0) > P(Qi(t) = Q 2 (t) = Q»(t) = 0) = ,g oo- The condition of ergodicity 
of the Markov chain (i.e., finite delay) is <?ooo > |2), from which ([T5T > follows. 

We now compute the end-to-end delay. Successively applying de l'Hopital's rule, and recalling ((28), the mean 
queue size at the source is found to be 



do 



xll 
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K" K" 
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(31) 



! 2^n 2K[ 

Moreover, recalling (|2H . the mean queue sizes at the first and second relays (or, equivalently, the busy probabilities 
since the size of the buffers is unity) are 
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(32) 



From © and CEO"©, (O follows. ■ 

In the particular case of a symmetrical system, i.e., pio,i = Pio,2 = Pio.3 = Pio an d pao.i = P20,2 = P20, the 
expression for the saturation throughput given in ([15} simplifies considerably. The result is stated in the following 
proposition. 



Proposition 3 The saturation throughput of a symmetrical three-hop system is 

T s = u(Pio,P20,Pii)v(pio,P20,Pnr 1 , 

where 

u(PlQ,P20,Pll) =P10(P 2 S P11 +P 2 S P10P11 +pIqPh) 



(33) 



v(piO,P20,Pll) =Pl0(l +P20)(PsPll +ploP20Pll) + (PlQ + P20) (PsPlOPll + P20Pll)- 
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Proof: Follows directly from (Q3} by setting p 10j i = pio,2 = Pw,a = Pw and p 2 o,i = P20,2 = P20- ■ 
Eq. ( f33b is amenable to interpretation for particular values of the parameters pio,P20,Pii- For example, letting 
P20 = 0, i.e., not allowing two-hop transmissions, yields 

Pw 



P11+P10P11 

For pn > 0, the denominator is < 3, which reflects the gain with respect to a system where intra-route spatial 
reuse is not permitted (pn = and r s = pio/3). 

The derivation of a closed-form expression for D from ( fT4t requires the constant C. We were not able to 
determine C analytically, but an approximation may be obtained as follows. Considering a symmetrical systemQ, 
and setting x\ — X2 = X3 = in d22b gives 

P10S001 +P2offoio = (-1 + A -1 +p w +P2o)gooo- 

Letting 5001 ~ Soirj. we have 

C *TT^Z T- (34) 

HPw +P20) 

Since g oi = P(Qi(i) = Qz{t) = 0) and g io = f(Qi(t) = Qz{t) = 0), we are approximating the probabilities that 
the first two queues are empty, and that the first and third queues are empty, as equal. Note that (|34| | is proportional 
to A. This is reasonable, since, for A — > 0, P(Qi(t) = Qa(t) = 0,Q3(t) = 1) — > or gooi — > <?ooo- The accuracy 
of the approximation is demonstrated with numerical results in the following section. 

IV. Numerical results 

We initially present numerical results for symmetrical two and three hop systems (Figs. [TJIUl and in Fig. [5] we 
examine a non-symmetrical three-hop system. The considered channel model consists of path-loss r~ a at distance 
r, where a is the path-loss exponent, and fading h which is constant within a slot, and spatially and temporally 
independent. We assume that h is exponentially distributed with mean one (i.e., \fh is Rayleigh distributed). The 
(instantaneous) received power is P r = Pr~ a h, where P is the transmit power, assumed common for all nodes. 
We define the average received SNR over a single hop as 7 = Pr~ a /No, where A*o is the thermal noise power. 
Assuming that a packet is successfully received if the received signal-to-interference-and-noise ratio (SINR) is larger 
than a threshold 6, the probabilities pw,P20,Pn defined in Section HT1 are 

pio = e- e '"<, p 20 = e- 2 " 9 ^, p n = p 10 /(l + 0). 

Apart from the considered "smart" opportunistic protocol (S-OPP), described in Section [Til for comparison 
purposes we consider the following two protocols. 

Multi-hop (MH): packets can only be transmitted over a single hop and nodes are divided in groups based on their 
spatial separation d = 1,...,N (in hops). In each slot, all nodes in a group can transmit simultaneously, and, 

'The approximation may be obtained easily for a non-symmetrical system as well. 
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across slots, a TDMA schedule is followed amongst the groups. If d = 1, all nodes can transmit in a given slot 
(full spatial reuse), while, if d — N, MH becomes a pure TDMA (round-robin) protocol. 

Regular opportunistic (OPP): The only difference between OPP and S-OPP is that if a node has a packet in its 
queue, it transmits, independently of the queue state of the successive node. 

Note that, in terms of feedback, MH only requires that the transmitter know whether its successive node 
successfully received the packet. In general, OPP and S-OPP require a more refined feedback, since a transmission 
has multiple potential receivers and all of them have to be informed of the outcome. This can be accomplished 
within a separate feedback slot, where, in a round-robin fashion, each node in the network (excluding the source) 
declares if it successfully received a packet and from which node. On the other hand, OPP and S-OPP do not 
require the scheduling of packet transmissions on which MH is based. 

For each protocol, we determine via simulation the average delay of the packets that arrive at the destination 
over a period of 10 6 slots. In the simulations, we relax assumptions (A1)-(A3), allowing for direct transmissions 
over distances exceeding two hops (if SINR > 6 is satisfied), taking into account interference from all transmitting 
nodes, and letting the relay buffer size B r > 1. The implication of B r > 1 is that a relay which has a packet in its 
buffer at time i, may receive a packet in slot t if it is silent. Unless otherwise stated, a = 3, 7 = 8 dB, 9 = 3 dB 
and B s — B r — 50, where B s denotes the source buffer size. 

In Figs. [T]|2j the delay is plotted vs. A for a two and a three hop system, respectively. Expectedly, S-OPP 
outperforms OPP and MH (for all possible d). Under little traffic, it is as aggressive as OPP, harnessing good fading 
conditions to perform direct two-hop transmissions. Under high traffic, it still behaves opportunistically, but avoids 
causing unnecessary interference, yielding a throughput gain of about 10% with respect to OPP under saturation. 
In fact, Fig. Q] depicts nicely how OPP suffers from interference for high traffic, resulting in larger delay than MH 
for A > 0.33. Note that the analytical approximation of the delay in Fig. |2]is satisfactory for all A. 

In Figs. [30 the simulated delay is plotted vs. A for four and five hop networks. The MH curves are obtained by 
selecting the delay optimal d for each A (which is d = 1 for the given system parameter values). For 7 = 5,10 dB, 
the maximum throughput of S-OPP is about 15% and 10% larger than OPP, respectively. For 7 = 5 dB in particular, 
the curves of OPP and MH are overlapping, due to the fact that two-hop transmissions are very rare. This implies that 
the gain of S-OPP with respect to OPP results only from the smart interference management. Another interesting 
observation is that the performance of S-OPP is insensitive to B r (seen by the light lines which correspond to 
B r = 1). The reason is that the events where three or more consecutive nodes have packets to transmit are quite 
rare for A smaller than the saturation throughput; therefore a relay buffer size larger than unity does not result in 
a notable end-to-end delay benefit. 

Closing the paper, we consider the performance of S-OPP in a three-hop system with non-equidistant relays. In 
Fig. |3J we plot the relay positions that maximize t s as a function of 7, and compare them with the respective ones 
obtained via simulation of a saturated system. Note that 7 in the non-symmetrical case is defined as 7 = 3 Q 7tot, 
where 7 to t is the end-to-end receive SNR. For normalization purposes, we set the source-destination distance to 
unity. If 7*1, T2 <E (0, 1), 7*1 < r-2, denote the distances of the first and second relays from the source, the probabilities 
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defined in Section [TT] are given by 

Pio.l = exp 
Pio.a = exp 

Pl0,3 = CXp 

Pu = 



(3ri) a 6» 
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and r s is obtained by 
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(35) 



(36) 



t s is evaluated as a function of 7 by substituting ( f35l > in d36*b . The discrepancy of the theoretical and simulated 
curves for 7 > 10 dB observed in Fig. [5] is due to the fact that, in the simulated system, the destination can be 
reached directly from the source with positive probability, which is not taken into account in the analysis. Focusing 
on the more realistic SNR range 6—10 dB, the main conclusion drawn from Fig. [5] is that, under normal S-OPP 
operation, it is advantageous to move the first relay slightly closer to the destination than 0.33. This position achieves 
the best tradeoff between reducing interference from the second relay and advancement towards the destination. 
This can be confirmed by the curves obtained when either two-hop transmissions or intra-route reuse are forbidden. 
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Figure 1, End-to-end delay versus A for a two-hop system. The solid line corresponds to D j5J and the dotted vertical line to 
that the simulation markers for S-OPP lie exactly on the theoretical curve, (a = 3, 7 = 8 dB, 6 = 3 dB, B s = B r = 50) 
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Figure 2. End-to-end delay versus A for a three-hop system. The solid line corresponds to the analytical approximation of D and the dotted 
vertical line to r s {33). (a = 3. 7 = 8 dB, 6 = 3 dB, B s = B r = 50) 
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Figure 3. End-to-end delay versus A for a four-hop system and 7 = 5, 10 dB. The light lines correspond to B r = 1 and the dark lines to 
B r = 50. (a = 3, 6 = 3 dB, B 3 = 50) 
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Figure 4. End-to-end delay versus A for a five-hop system and 7 = 5, 10 dB. The light lines correspond to B r = 1 and the dark lines to 
B r = 50. (a = 3, 6 = 3 dB, B 3 = 50) 




Figure 5. Throughput-optimal relay positions (normalized to unity) for a three-hop S-OPP system vs. equivalent link-SNR 7. The dark lines 
are obtained by maximizing the theoretical saturation throughput and the light lines by simulation, (a = 3, 9 = 3 dB, B s = B r = 50) 



